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We present here the latest results from the QCDSF collaboration for moments of gener- 
alized parton distributions and transversity in two-flavour QCD, including a preliminary 
analysis of the pion mass dependence. 



1. INTRODUCTION 

For many years deep-inelastic scattering experiments have provided a wealth of infor- 
mation regarding the quark and gluon content of the nucleon, mainly through parton 
distribution functions which describe the longitudinal momentum distributions of quarks 
and gluons in the nucleon. Although the amount of knowledge gained from such experi- 
ments is constantly increasing, little is known about the transverse structure and angular 
momentum distribution of partons within the nucleon. Generalized parton distributions 
(GPDs) [ H] have opened new ways of studying the complex interplay of longitudinal 
momentum and transverse coordinate space [El Ej, as well as spin and orbital angular 
momentum degrees of freedom in the nucleon [ Ej . A full mapping of the parameter space 
spanned by GPDs is an extremely extensive task which most probably needs support from 
non-perturbative techniques like lattice simulations. 



*Talks presented by Ph. Hagler and J. Zanotti at Baryons 2004. 
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Besides the polarized and unpolarized GPDs, there exist also four independent ten- 
sor/helicity flip GPDs HT{x,^,t), Et, Ht and Et, as has been shown in The first 
one, Hxixj^jt), is called generalized transversity, because it reproduces the transversity 
distribution in the foward limit, Ht{x, 0,0) = 5q{x) = hi{x). Since the quark tensor 
GPDs flip the helicity of the quarks, they do not contribute to the deeply virtual Comp- 
ton scattering (DVCS) process 7*p •yp'. Although this could in principle be balanced 
by the production of a transversely polarized vector meson, 7*p —>■ uitp' ■, it has been 
shown that the corresponding amplitude vanishes at leading twist to all orders in pertur- 
bation theory [ 13 E] • The only process we are currently aware of which gives access 
to the generalized transversity is the diffractive double meson production proposed in [ 
IH]. It will be very interesting to see if a measurement of tensor GPDs in this process is 
indeed feasible. In any case, the situation seems to be much more difficult as compared to 
the (un-) polarized GPDs, and an independent lattice calculation of the moments of the 
helicity flip GPDs would be highly valuable. 

There has been a large amount of activity within the lattice community in the area 
of GPDs leading to some exciting results from the QCDSF [ HUJ HU IH 113 and LHP [ 
d 113 CHI im im Hi collaborations. 

Here we present the latest results from the QCDSF collaboration for the first three 
moments of the unpolarized and polarized GPDs and results for the tranversity GPDs. 
All results are preliminary. 

2. SIMULATION DETAILS 

We simulate with N{ = 2 dynamical configurations generated with Wilson glue and non- 
perturbatively 0{a) improved Wilson fermions. For five different values jS = 5.20, 5.25, 
5.26, 5.29, 5.40 and up to three different kappa values per beta we have in collaboration 
with UKQCD generated C(2000 — 8000) trajectories. Lattice spacings and spatial volumes 
vary between 0.075-0.123 fm and (1.5-2.2 fm)'^ respectively. 

Correlation functions are calculated on configurations taken at a distance of 5-10 tra- 
jectories using 8-4 different locations of the fermion source. We use binning to obtain 
an effective distance of 20 trajectories. The size of the bins has little effect on the error, 
which indicates auto-correlations are small. This work improves on previous calculations 
by adding one more sink momentum, p2, and polarization, Fi. We use po = (0,0,0), 
Pi = (p,0,0), p2 = (0,p,0) {p = 27r/Ls) and Fu^poi = ^(1 + 74), Fi = i(l + 74)^7571, 

r2 = 1(1 + 74)^7572- 

3. GENERALIZED PARTON DISTRIBUTIONS 

For a lattice calculation of GPDs, we work in Mellin-space to relate matrix elements 
of local operators to Mellin moments of the GPDs. For quark distributions, the twist-2 
operators are 







= qi^^' D ■■■ D q , 


(1) 






= qi^T^^ D ■■■ D q , 


(2) 




■■/in} 


= gz7V^^i7^2 D ■■■ D g , 


(3) 
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Figure 1. Generalized form factors A20, -B20 ond C2, together with a dipole fit for the 
non-singlet, u — d (left), and singlet, u + d (right). 

where D= |(-D — D) and {■■■} indicates symmetrization of indices and removal of 
traces. The non-forward matrix elements of the twist-2 operators Eqs. |21 El) specify 
the {n — ly^ moments of the spin-averaged, spin-dependent and spin-flip generalized 
parton distributions, respectively. In particular, for the unpolarized GPDs, we have 

j'^dxx^'-'H'^{x,^,A') = Hlii,^^) , 

1' rfxx"-^E^(x,e,A2) = Elii.l^^) , (4) 
where [H] 

n-l 

i^^(e,A^) = ^<2,(A')(-2ef + Mod(n+l,2)C^(A^)(-2er, 

1=0 

n-l 

El{i.A') = X]<2.(A^)(-20''-Mod(n + 1,2)C^(A2)(-20", (5) 

and the generalized form factors A^2j(A^), B'^2ii'^'^) ^n(A^) for the lowest three 
moments are extracted from the nucleon matrix elements {p'\0^'^^"'^"^p) [H]. Note that 
the momentum transfer is given by A = p' — p with t = A^, while ^ = —n ■A/2 denotes 
the longitudinal momentum transfer. For the lowest moment, Aio and Bio are just the 
Dirac and Pauli form factors Fi and F2, respectively. We also observe that in the forward 
limit (A^ = ^ = 0), the moments of Hg reduce to the moments of the unpolarized 
parton distribution A^o = (x""^). Finally, the forward limit of the x-moment of the GPD 
E, J dx X E{x,0,0) = i?2o(0), allows for the determination of the quark orbital angular 
momentum contribution to the nucleon spin, L'^ = 1/2 ((x) -|- B20 — Ag), where (x) is the 
quark momentum fraction 

In order to extract the non-forward matrix elements we compute ratios of three- and 
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Figure 2. Generalized form factors A 
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together with a dipole jit. All form factors have been normalized to unity. 
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(6) 



up to known kinematical factors as long as <^ r <^ t < \^t- 

In Fig. n we show, as an example, the generalized form factors A20, -B20 and C2 in the 
MS scheme at /i^ = 4 GeV^ for the non-singlet, u — d (left), and singlet, u + d (right), 
on a 24^ x 48 lattice at (3 = 5.40 and Hgea = ^^vai = 0.13500 corresponding to a lattice 
spacing, arg = 6.088 and = mpg ^ 970 MeV. We note here that in the calculation 
of the singlet matrix elements, we neglect contributions coming from disconnected quark 
diagrams as these are extremely computationally demanding. 

The generalized form factors A^q"^, ^20'^ -820"'^ ^^^^ described by the dipole 
ansatz 

while C2 ""^ is consistent with zero. Our result for Bl^'^'^ reveals a small but positive value 
at this heavy quark mass. Future work will include a full non-perturbative analysis of 
^qu+d Ag"~°', together with a chiral extrapolation, to reveal the quark contributions 
to the nucleon's total spin and orbital angular momentum. 

Burkardt [0] has shown that the spin- independent and spin-dependent generalized par- 
ton distributions H{x, 0, A^) and H{x, 0, A^) gain a physical interpretation when Fourier 
transformed to impact parameter space at longitudinal momentum transfer ^ = 

= / 1^ e-^^--^-H{x, 0, -Ai) , (8) 
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Figure 3. Generalized form factors A^q"^, ^20''^ ^''^^ ^scT"'? together with a dipole fit. All 
form factors have been normalized to unity. Dipole masses are extrapolated linearly to the 
chiral limit. 



(and similar for the polarized Aq{x,b±)) where q{x,b±) is the probability density for 
a quark with longitudinal momentum fraction x and at transverse position (or impact 
parameter) b±. 

Burkardt [ 3j also argued that H{x, 0, — A^^) becomes A^-independent as x ^ 1 since, 
physically, we expect the transverse size of the nucleon to decrease as x increases, i.e. 
lim^;^! q{x, b±) oc S^{b±). As a result, we expect the slopes of the moments of H{x, 0, — A^) 
in A^ to decrease as we proceed to higher moments. This is also true for the polarized 
moments of H{x, 0, — A^), so from Eq. (0) with ^ = 0, we expect that the slopes of the 
generalized form factors yl„o(A^) and yl„o(A^) should decrease with increasing n. 

In Fig. 121 we show the A^-dependence of Ano{A'^) (left) and A„o(A2) (right), n = 1, 2, 3, 
for f3 = 5.40, Ksea = i^vai = 0.13500. The form factors have been normalized to unity 
to make a comparison of the slopes easier and as in Fig. ^ we fit the form factors with 
a dipole form as in Eq. ((7j). We observe here that the form factors for the unpolarized 
moments are well separated and that their slopes do indeed decrease with increasing n 
as predicted. For the polarized moments, we observe a similar scenario, however here 
the change in slope between the form factors is not as large. The flattening of the GFFs 
Ano(A^) has first been observed in Ref. [El, where at the same time practically no change 
in slope has been seen going from A2o(A^) to A^o^A'^). 
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Figure 4. The iso-vector tensor charge vs. m^. 
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Although fitting the form factors with a dipole is purely phenomenological, (see Ref. [E 
for an alternative ansatz) it does provide us with a useful means to measure the change in 
slope of the form factors by monitoring the extracted dipole masses (Mi, M2, M3) as we 
proceed to higher moments. We have calculated these generalized form factors on a subset 
of our full complement of (/3, k) combinations and have extracted the corresponding dipole 
masses. We plot these dipole masses in the lower-right of Fig.Elas a function of rrip^. The 
values quoted in the chiral limit are calculated via a combined dipole-mass/pion-mass fit 



^dipole, 



(9) 



depending on the two parameters and a, in order to extract the maximum amount of 
information from our available (/3, k) combinations. Having done the fit, we may shift the 
raw numbers to a common curve given by Eq.ijH)) at the physical pion mass, m7r,phys = 139 
MeV. The results of this procedure are shown in Fig. El for the GFFs Ano, n 



1,2,3. 



The important feature to note in Fig. El is the distinct separation between (and increase 
in magnitude of) the dipole masses as we move to higher moments {x —>■ 1). Although the 
data available for M3 is limited, the behaviour of all dipole masses appears to be linear 
with m|g. Consequently, we perform individual linear extrapolations of the dipole masses 



Ml, M2, M3 to the physical pion mass, although the findings of Ref. [|2l| suggest that 
the chiral extrapolation of the dipole masses of the electromagnetic form factors may be 
non-linear. 

If the dipole behavior Eq. ((7|) continues to hold for the higher moments as well, and if 
we assume that the dipole masses continue to grow in a Regge-like fashion, we may write 



with Mf = const. + l/a', where const. ^ —0.5 GeV^ and 1/a' ^ 1.1 GeV^. This is 
sufficient to compute Hq{x, 0, A^) by means of an inverse Mellin transform [I22j. 

Having done so, the desired probability distribution of finding a parton of momentum 
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Figure 5. Combined dipole-Zpion-mass fits to the generalized form factors Atio C'l^d At2o- 
fraction x at the impact parameter b± can then be obtained by the Fourier transform of 



We now turn our attention to our latest results for the first two moments of the tensor 
GPDs. Since the parameterizations for the Mellin-transformed tensor GPDs have been 
derived recently in [|2Sl I21| in terms of the GFFs Axni, Axni, Btui and -Brm) we are now 
in a position which allows for the extraction of the moments of Hx{x,^,t) etc. from our 
lattice calculation. The lowest moment of the generalized transversity, J dxHT{x,^,t) = 
^^10(^^)5 is equal to the tensor form factor (7t(A^) = AtioIA^). Our results for the iso- 
vector tensor charge, g^i^), are shown in Fig. |3]as a function of the pion mass squared. 
The lattice results for g^ have been non-perturbatively renormalized and transformed to 
the MS scheme at = 4 GeV^. The plot indicates that sizeable discretization and/or 
finite volume effects may be present which we will study in detail in the near future. 
A linear extrapolation of the isovector tensor charge to the physical pion mass gives 
5't(0) — 1-09 ^ 0.02. A more advanced form for the chiral extrapolation of g^ including 
effects from the pion cloud is also currently under investigation [ 12^1 12n] • 

Fig. shows combined dipole-/pion-mass-fits of (yf|.(A^) and the n = 2 -moment of the 
generahzed transversity, At2o[^^), following Eq.(jni). The corresponding dipole masses 
and fit-parameters are given by 

Atio ■ M° = 1.63±.02 GeV, a = 0.10±.03 GeV"^ xVDOF = 1.1 

At2o ■■ = 2.17±.04 GeV, a = -0.06±.06 GeV"^ xVDOF = 0.9, (11) 

showing an increase in the dipole mass, going from the n = 1 to the n = 2-Mellin-moment. 
The errors in Eq. (fTT|) are purely statistical. Results for the other tensor GFFs Atuo, Bxno 
and Bttio for n = 1,2 will be presented in a future publication. For the tensor GPDs, the 
number of GFFs grows rapidly with ra, i.e. A^ = 2\n/2 \ + n. For = 3 we would have 
to extract simultaneously seven independent GFFs from the lattice, which poses quite a 
challenge. We plan to study the relevant two- derivative tensor operators on the lattice to 
see if and to what extent an extraction is possible. 



Eq. dHD. 
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